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Abstract 



We prove that the periodic initial value problem for a modified 
Euler-Poisson equation is well-posed for initial data in ff*(T™) when 
s > m/2 + 1. We also study the analytic regularity of this problem and 
prove a Cauchy-Kowalevski type theorem. After presenting a formal 
' derivation of the equation on the semidirect product space Diff k C°° (T) 

^ , as a Hamiltonian equation, we concentrate to one space dimension 

' (m = 1) and show that the equation is bihamiltonian. 

In this paper we study the periodic Cauchy problem for the modified 
, Euler-Poisson equation (mEP) 

in 

^ ! dtn + div(ni;) = , . 

^ ; dtv + {v.V)v + gradA-^n = 0, x G T™, f G M. ^™ ' 

^ I as well as its hamiltonian structure and integrability. 

^ ' The equation (mEP) is related to the Euler-Poisson equation 

^ ■ dtn + div(nv) = 

dfV + {v.V)v + gradcj) = (1) 
A(p - e't' + n = 0,x eT"',t £R 

which describes the fluctuations in the ion density of a two-component 
plasma of positively charged ions and negatively charged electrons (therefore 
it is also called ion acoustic plasma equation |LiSatj ^. 

*to appear in Transactions of the American Mathematical Society (in press) 
°2000 Mathematics Subject Classification: 35Q53, 35Q05, 35A10, 37K65. 
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Linearizing the operator A^((/>) = e"^ — /S.cj) at (/> = in the Euler-Poisson 
equation , we obtain the local form of the modified Euler-Poisson equation 
(mEP): 

dtn + div(ni;) = 0, 

dtv + {v.V)v + grad(/) = 0, (2) 
A0-(/. + n = O, xGT"", t GM. 

The equation @ , like the Euler-Poisson equation , admits an approxima- 
tion which preserves dispersion and leads to KdV (see remark Inverting 
the operator A := / — A, we write the system in Q in the nonlocal form 
(mEP). 

Besides its relation to the Euler-Poisson equation , the modified Euler- 
Poisson equation is also remarkable for its bihamiltonian structure in one 
space dimension that we describe here. 

Many different techniques have been developed based on Picard's con- 
traction theorem on Banach spaces in the study of nonlinear partial differ- 
ential equations. One approach originated in an observation of V. Arnold 
|Arnj that the initial value problem for the classical Euler equations of a 
perfect fluid can be stated as a problem of finding geodesies on the group 
of volume preserving diffeomorphisms. Subsequently, this observation was 
used by D.G. Ebin and J. Marsden in |EMa| who developed the necessary 
functional analytic tools and established sharp local well-posedness results 
for the Euler equations in a class of Sobolev spaces. 

The first section of this work is devoted to develop an appropriate an- 
alytic framework for the modified Euler-Poisson equation (mEP) using a 
similar approach and prove the following theorem. 

Theorem 1. For s > ni/2+1, given any initial data (no,Wo) ^ //^""^(T™, M) x 
iJ'*(T™, M*"), there exists aT > and a unique solution {n,v) to the Cauchy 
problem for the modified Euler-Poisson equation (mEP) such that 

V G C{[0,T),H'{T"',R"'))nC\[0,T),H'-\T"',R"')) 

and 

n G c([o,r),i/"~i(T™,M)) ncH[o,r),i/"-2(T",M)) 

and the solution {n,v) depends continuously on the initial data (no,fo)- 

Another powerful tool in the study of partial differential equations is 
the Cauchy-Kowalevski theorem. An abstract version of this theorem was 
developed by L.V. Ovsjannikov |Uvsl[ IDvs2j . F. Treves |Trej . L. Nirenberg 
|Nirj . T. Nishida |Nisj and M.S. Baouendi and C. Goulaouic |BGj among 
others and subsequently applied to the Euler and Navier-Stokes equations. 
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The study of analytic regularity of solutions to the Camassa-Holm equa- 
tion by A. Himonas and G. Misiolek |HMlj . |HM3j using this abstract theo- 
rem led us to investigate the analytic regularity for (mEP). In Section 2, we 
prove the existence and uniqueness of local analytic solutions to the Cauchy 
problem for the equation (mEP). 

Theorem 2. // the initial data (no,fo) is analytic on T™ x T™ then there 
exists an e > and a unique solution (n, v) of the Cauchy problem for the 
equation (mEP) that is analytic both in x and t on T™ x for all t in 



This result can be viewed as a Cauchy-Kowalevski type result for the 
equation (mEP). Even though the equation (mEP) admits an approxima- 
tion by the Korteweg-De Vries equation, the analytic regularity results for 
the two equations are quite different. In contrast with the Korteweg-De 
Vries equation whose solutions are analytic in the space variable for all time 
but not analytic in the time variable (see |Truj . |KaMj ) . the solutions to 
the modified Euler-Poisson equation are analytic in both space and time 
variables. 

In the third section we derive the equation (mEP) as a Hamiltonian 
equation on the semidirect product space Diff(T™) tx C°°(T™) following the 
treatment of V. Arnold and B. Khesin in |AKj and J. Marsden, T. Ratiu 
and A. Weinstein in ,MIlWi of the Hamiltonian formalism related to fluid 
and gas dynamics. Then we concentrate to the one space dimension m = 1 
and prove the following theorem. 

Theorem 3. For m = 1 the modified Euler-Poisson equation (mEP) is 
bihamiltonian with the pair of hamiltonian functionals 



In the proof we use prolongations to check the compatibility of the in- 
duced Poisson brackets by these hamiltonian structures. In particular the 
modified Euler-Poisson equation (mEP) can be derived as a Hamiltonian 
equation on the semidirect product space of the Virasoro algebra with the 
smooth functions on the torus vir x C°°(T) along with a nonlocal hierarchy 
of equations called Hunter-Zheng equations (see |BDPj for the bihamiltonian 
structure of the Hunter-Zheng equations). 




and 




nv dx. 



3 



Remark 1. The Korteweg-de Vries equation (KdV) can be derived as an 
approximation to the Euler-Poisson equation by a perturbation analysis (see 
\Sa^ ). Using this approach it is straightforward to obtain an approximation 
to the system of equations in ^ which preserves the dispersion and leads to 
KdV. 

1 Local well-posedness in Sobolev spaces 

In this section we study the Cauchy problem for the modified Euler-Poisson 
equation (mEP) where n = n{t, x) : M x ^ M and = v{t, j;) : M x ^ 
and A"^ = (/ _ A)-^ is the Bessel potential. 

In order to prove Theorem ^ we use the method of first restating the 
problem as an initial value problem for an ordinary differential equation on 
the group of diffeomorphisms of Sobolev class and then applying the 
existence theorem for vector fields on Banach manifolds. 

Proposition 1. For s > m/2+1, apair{n,v) G i?*-i(T",M)xi^'*(T™,M™) 
is a solution to the Cauchy problem for (mEP) with initial data {nQ,VQ) if 
and only if v = r] o j^'^ and n = C,o 7^-*^ where (7, C,, rj) is a solution to 

dtC = -((C o l~^)dw{r] o 7-I)) o 7, 

at77 = -(grad(A-2(Co7-i)))o7, (3) 

da = V, 

with initial conditions C(0,x) = nQ{x),r]{0, x) = vq{x),^{0, x) = idx- 

Therefore the Cauchy problem for (mEP) can be reformulated as an 
initial value problem for the ordinary differential equation 

|(7,C,^) = (r/,F(7,r/, 0,^(7,0) (4) 

where 

G(7,C) = -(gradA-2(Co7-i))o7, 
C, V) = -((C o 7"^)div(7? o 7-I)) o 7) 

with initial data (70,%, Co) = {'i'dx,vo{x),no{x)). 

In the proof of Theorem ^ we repeatedly use three standard results on 
Sobolev spaces: The Schauder ring property, Sobolev imbedding theorem 
(that we refer to as Sobolev lemma) and the composition lemma (see, for 
example, [Sj and |Adj ) . 
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Proof of Theorem ^ If the map 

ViT"") X ^r'*-i(T'",M) X i?'^(T™,M™) iJ*(T'",M'") x H'-'^(J"',R) x i?'^(T",M'") 

is locally Lipscliitz then by the fundamental theorem for ordinary differential 
equations on Banach spaces |Dij there is a unique solution 

(7,r/,C) G ViJ"") X H'iT"',R"') x H'-\T"',R) 

to the problem Q for s > m/2 + 1 with initial data 

C(0,x) =no(x), ??(0,x) =t;o(x), 7(0,3;) = x. 

Note that the dependence of the solution of problem Q on initial data 
is smooth. However the map 7 1— > 7^^ on is continuous but not C^. 
Therefore we only have continuous dependence on initial data of the solution 
to the Cauchy problem for (mEP). 

By Proposition n the proof of Theorem ^ is reduced to showing that the 
maps 

7 ^F(7,r/,C) G L{H%T"',R"'),H'-^{T"\R)), 
C ^F(7,r/,C) e L{H'-\T"',R),H'-\T"',R)), 
r]^F{-f,r]X) £ L{H'{T'^,R"'),H'-\T"'\R)), (5) 
7 ^G(7,C) G L{H'iT"',R'^),H%T"',R"')), 
C ^G(7,C) G L{H'-\r"'\R),H'{T'',R^)) 

are locally Lipschitz in 7, and C (uniformly with respect to the remaining 
variables). 

In the following estimates, the subscripts 7 and C of a constant indicate 
the dependence of the constant on ||7||hs and ||C||h''-i respectively. 

7 I— > G(7, C) is locally Lipschitz: 

Let 71,72 G V^T"") and C G H'~'^{T'^ ,R). By the composition lemma, 
||gradA"2(Co7fi) 071 - gradA~^(C o 7^^) o-/2\\h'' 
< a,, ||gradA-2(^ o 7-I) 07- gradA-2(^ o j^^)\\h^ 

where 7 = 71 o 72"^- Let C = C ° ■ Then it is enough to show that the 
following estimate holds: 

||gradA"2(Co7"^) o7-gradA"2CI|//= < C^,(;h - idxWn- ■ (6) 
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Next we show that holds for s > m/2 + 1. The left side of the 
inequality © has the following form 



llgradA 2(^0 7 i)o7-gradA ^^H^, 

~ ||gradA"^(C o 7"^) 07- gradA~^CIlL2 

+ ||div(gradA-2(^ o 7-I) 07- gradA-2^)||^.-i. 



(7) 
(8) 



We first estimate the term in ((71). For any r > 111/2 + a we have the 
Sobolev imbedding into a Holder space 

H'-iT"') ^ C7'^(T™) 

with a bound 

\u{x) - u{y)\ < C\\u\\Hr\x - yl" , (9) 

for any x,y G T. Therefore using the composition lemma and applying Q 
with r = s — 1 and a = {s — 1 — m/2)/2 we find 

|[gradA~^(C o 7"-^) 07 — gradA^^CllL2 

< ||gradA~^(C o 7"^) 07- gradA"^(C o 7~"^)||l2 + ||gradA"^(C o 7"^) - gradA'^ClliZ 

< C^\\C\\h^-i (IIt - ^411?^' + ||7~' - id^^Hs) . 



Adding and subtracting the appropriate terms we estimate the ^ term 
in (jH)) by the sum 



< 



5^5^a,-5,A-2(Co7-i)o7(9oi-5D 
i=i j=i 

+ \\AA-\C o 7'^) ° 7 - AA-2(||^.-i. 



(10) 



(11) 



Using Schauder ring property and composition lemma the first summand in 
(jlOp is bounded by 



< IICII//=-i||7 - id.. 



In order to estimate the second summand (jllj) we add and subtract the 
terms A^^((^ o 7^"*^) o 7 and A^^(^. After cancellations we obtain 



= ||A-2(Co7-i)o7-A-2c||^.-i 
<||A-2(Co7-i)o7-A-2(Co7-i)||^._i 

+ ||A-2(Co7-i)-A-2Ck^-i- 



(12) 
(13) 
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Let u be A'^^C Then we have 

||A-2(Co7-i)o7-A-2(Co7-i)|[^.-x = 

< Ikllci IIt - idxWn''-^- 

Using the Sobolev lemma with the composition lemma we obtain the esti- 
mate 

<C7llCll//-i|l7-^4||H^ (14) 

for (Ha). 

For s > m/2 + 2 the term H13() can easily be estimated like (|12() . For 
m/2 + 1 < s < m/2 + 2 we first observe that the estimate 

\\A'\C o 7"') - a-'CIIh-i < lie o - CIIl^ (15) 

holds for s — 3 < and then using Q as before, we obtain the following 
estimate for ((Ti^ 

<C4C\\H^-4T'-^dxrHs (16) 

where a is equal to (s — 1 — m/2)/2 > 0. However the assumption s — 3 < 
does not follow from s < m/2 + 2 if m > 3. Nevertheless one can use the 
following inductive argument until s — {2k + 1) < (it ends in finitely many 
steps since s < m/2 + 2). 

If s — 3 > we split (|13|) as in 0-®. The L2 part can be estimated as 
for s - 3 < 0. The H^"^ part 

= ||gradA"2(C o 7"^) o 7.1)7 - gradA-2c||iy.-2 (17) 

is bounded by 

< ||gradA-2(C o 7-I) o 7.P7 - 1) \\hs-2 (18) 

+ |lgradA-^(C o 7"^) 07- gradA-2^||^^.-2 (19) 

Here the first summand (|18|) is estimated using the Schauder ring property 
with the composition lemma 

||gradA-2(Co7-i) 07.(^7- 1)||^,_2 < (7^||C||//s-i II7 - i^lli/- 



/ Du{y)dy 

J X 
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For the second summand (fT^ we use the steps (|H])-(fTl]) to reduce it to 
estimating 

\\A-\Co7-')-A-'C\\h^-^. (20) 

If s — 5 < we proceed as in ()15|) - ()16() . Otherwise we repeat the steps 
(fT7|)-(PH). 

7 1-^ F{j, C, rj) is locally Lipschitz: 

Let 71,72 G i7"(T™,]R'^) and C G i^'^^T™, M), r/ G H%T"',M."'). Then 
by Schauder ring property we have 

||C(div(77 o 7-I) o 71) - C(div(?7 o 7-^) o 72)||h''-i 

< \\C\\H''~A\^dmVjdj{ji^)m0 1l - dmrijdj(j2^)m°l2\\H'-^ 
^ C^i,r2,(^\\9mVj\\H--A\dji7i^)m - dj{j2^)m\\H--i. 

Using the Schauder ring property one more time we bound this term by 

^ C'^,i,12,C,v\\v\\h^\\1i -12\\h=- 

and therefore 7 1— > -F(7, C; v) is locally Lipschitz. 

It is straightforward to show that the second, third and fifth maps in © 
are uniformly Lipschitz using properties of Sobolev spaces. This completes 
the proof of Theorem ^ 

□ 

Next we observe that the Cauchy problems for the equation (mEP) and 
the Euler equations of an incompressible fluid are not only similar for low 
regularity (Sobolev class) data but also for high regularity (analytic) data. 

2 Analytic regularity 

In this section we give a proof of theorem|2lthat states the analytic regularity 
(i.e., existence and uniqueness of analytic solutions for analytic initial data) 
of the Cauchy problem for (mEP). 

Our approach is motivated by the work of M.S. Baouendi and C. Goulaouic 
|BGj who studied analytic regularity of the Cauchy problem for Euler equa- 
tions of incompressible fluids. 

The proof of theorem |^ relies on a contraction argument in a decreasing 
scale of Banach spaces Xg (i.e. if s' < s implies Xg C Xg' and ||| • \\\s' < 



For s > 0, let the spaces Eg be defined as 
E, = \u£ C^^iJ"") : f udx = Oand \\\u\\\s = sup ;7f;7,","; < oo 



||9^n||i^<Tsl'''l 



where a is any integer such that a > 1 + m/2 and let Xg be given by the 
Cartesian product Eg x Eg. The norm ||| • \ \\xs can be chosen to be any of 
the standard product norms on Eg x Eg. The following lemma states the 
ring property for the spaces Eg. 

Lemma 1. Let < s < 1. There is a constant c > which is independent 
of s such that we have 

lll^it^llls < Ills 

for any u,v G Eg. 

We omit the proof of lemma ^ (see |HMlj for the case m = 1). First 
we rewrite the equation (mEP) in a more convenient form. Let n and v be 
denoted by ui and U2 respectively. Then we can write the equation (mEP) 
in terms of (ui,U2) as 

dtui = Fi{ui,U2) := V2{uiU2), 

dtU2 = F2{ui,U2) := Vi{u2)u2 + TiVsUi 

where 

Viin) := — grad(n), ^3(71.) := A~^n, 

P2(u) := -div(t;), 'Pilu)v := -(V„n) = -{Du)v. 

The following lemmas give the suitable bounds on these operators to 
prove Theorem 121 

Lemma 2. For < s' < s < 1, we have 

C 

|||7'in|||s' < tIII^-IIIs, 

s — s' 

1117'2'yllls' < — ^illwlllc. 



Proof. By the definition of |||-|||s, we have 

^d^Vin\\H^s'\''\ 



|||Pin|||s' = sup 



N>o k\/{\k\ + l) 



9 



The norm on the right hand side can be written in the local coordinates 
up to a constant as 

m 

||a''(grad(n))||j^<T(Tr-,R-) ^ ^ \\d''djn\\H'y{T"^fi,). 

Then we have the estimates 

Pin < c sup sup 

|fc|>o|/3|=i k\/{\k\ + lY 

\\d^+Pn\\H-s\^\+^ {k + (5)\ s'\^\ {\k\ + lf 



' S {k + P)\/{\k\ + 2)2 {\k\ + 2)2 .Ifcl+i A;! • 
Note that 

Therefore we have 

WmWl, < m n U sup sup |,|,/ \^ . 22 



Note also that 



Then it follows from the formula ((23} that the inequality 



(k + BV 

sup ^ = sup (A:^ + 1) < \k\ + 1. 

|/3|=1 fe' l<i<m 



s'l'^l (fc + /3)! /|A;| + lV 



holds. By ((121) and we obtain 

C 

|||7'in|||s/ < -|||n|| 

s — s 

The estimate for 7^2 follows similarly. 

Lemma 3. For any < s < 1, the estimate 

\\\r3{u)\\\s < \\\u\\\s 

holds if u E Eg. 



< —7 (23) 



□ 
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Lemma 4. For < 5' < s < 1, we have 

|||^4('u)^^|||s' < rlll^'llls'ill'"! 



Proof. We write 'P4{u)v in terms of the linear operator Du as 'P4,{u)v = 
{Du)v. Then by Lemmas we have 

\\\Vi{u)v\\\s' = \\\{Du)v\\\s' 

m 

= III VidiUj\\\s' 

m 

< lll^i|||s'|l|5i%llls'- 

Now we reduce the proof to the case that we handled in the proof of lemma 

HI 

m 

|||^4(^i)t'|||s' < clllullls' ^ IPilijllls' 

^ \\d'+%\\H^s'\>'\ 

< c\\\v\\\s' / ^ sup sup 



,ttl?>o 1^1=1 W(l^l + 1)^ 

/ s'l^^i {k + py. f\k\ + i 

< C \ V s' U s sup sup -TT-T I, I , „ 

|fc|>o|/3|=i k\ \\k\ + 2 
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Clearly, to finish the proof, it is enough to show that 

(A: + l)<^. (24) 



5*^+1 \k + 2^ 

Let s' = As, < A < 1 and /(A) = (/c + 1)(1 - A)A^. Then, 



:/(A) 



s(l - A)' ' ' \k + 2j 

For A; = it is clear that /(A) < 1. For k > 1 the function /(A) = 
(A; + 1)(1 — A)A'^ is continuous in the interval < A < 1 and it has zeros 
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at the endpoints of the interval [0, 1] and a maximum at A = -^r^ such that 
/(fcTi) = (kTi) < 1- ^^^"^ ^^^^ 



k + 2 



{k + l)< 1. 



Therefore 



{k + l)< 



gk+i \k + 2j ^ ^ - s{l-X) 

and the formula H24|) holds. Using the formula 1)24^ we obtain the desired 
estimate 

|||^4('u)^^|||s' < rlH'wIlls'llkllls- 

s — s' 

□ 

Now we are ready to prove theorem |21 

proof of theorem |21 We refer to the version of the abstract Cauchy- 
Kowalevski theorem in ;Nis] . We only need to verify the first two conditions 
of this theorem since the map F{ui,U2) does not depend on t explicitly. 

Clearly, t i — > F{t, u{t)) = {Fi{ui, U2), F2{ui, U2)) is holomorphic if t 1 — > 
ui{t) and t I — > U2{t) are both holomorphic. We only need to show that 
Fi{ui,U2) and ^2(^*1, 1*2) are in Eg' if ui,U2 G Eg. We begin with estimates 
on Fi. 

By Lemma 121 and Lemma we have 

|||Fi(ni,n2)|||s' = |||'P2('"l'"2)|||s' < 7llhl|IW||li2|||s. 

s — s' 

Similarly, for F2 , using the lemmas |21 El and jll we have 

\\\F2{ui,U2)\\\s' = \\\r^{u2)u2+VlV3Ui\\\s' < ( 1 1 1^2 1 1 Is + 1 1 1 1 Is) ■ 

We proceed to establish the second condition of the abstract Cauchy- 
Kowalevski theorem. We will show that for some c independent of t, 

|||Fi(ni,U2) - ^1(^1,1^2)1118' < i\\\u - v\\\x, 

s — s' 

and similarly 

|||F2(ni,M2) - ^^2(^^1,^^2)1115' < :\\\u - v\\\x, 

s — s' 

where u = (ui,U2) and v = (^1,^2). 
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To obtain the first estimate above, after applying Lemma |21 we add and 
subtract the term U1V2 and use Lemma ^ 

|||Fi(ni,U2) - Fi{vi,V2)\\\s' = \\\'P2{uiU2 - ViV2)\\\s' 

c 

< 7 ( pi s Ki2 - ^^2 s + Fi ~ ^1 U K'a 

s — s' 

Then, assuming that |||ii|||s < R and |||w|||s < R, we have 

\\\Fi{ui,U2) - Fi{vi,V2)\\\s' < j\\\u - v\\\xs- 

s — s 

To estimate the F2 component, we use lemmas |21 and 01 

\\\F2iui,U2) - F2{vi,V2)\\\s' = \\\i'P4,U2)u2 - {ViV2)v2 + VlV:i{ui 

Q 

< \\\(.'P4U2)U2 - (V4V2)v2\\\s' H tH^I " 

s — s 

(25) 

Note that 

V4:{U2)U2 - 'P4{V2)V2 = V4{u2){u2 - V2) + 'Pa{u2 - V2)v2. 

Using Lemma 0] and the above identity, ()25() implies 

Q 

|||F2(ni,U2) - F2{vi,V2)\\\s' < 7 (|||n2|||s|||li2 - ■V2|||s' + 

s — s' 

+ |||'«2 - 'y2||U||w2|||s' + llllil - Vi\\\s) . 

Therefore the estimate 

\\\F2{ui,Vi) - F2{u2,V2)\\\s' < J^\\\U - v\\\x. 



holds. This completes the proof of theorem |21 



□ 



3 Bihamiltonian structure and integrability 

A number of partial differential equations that describe fluid motion can 
be derived as equations for geodesies on various infinite dimensional Lie 
groups. For instance, the Euler equation for ideal incompressible fluid flow 
is the geodesic equation on the group of volume-preserving diffeomorphisms 
of a Riemannian manifold M with a right invariant metric given by the 
inner product on the tangent space at the identity of the group |EMaj . 
Other examples are 
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• Korteweg-de Vries equation and Camassa-Holm equation on the Bott- 
Virasoro group (see for example |UKj and |Mis3j ). 

• ideal incompressible MHD (magnetohydrodynamics) on the semidirect 
product of volume preserving difFeomorphisms with the divergence free 
vector fields, 

• Hunter-Saxton equation on the homogeneous space of all diffeomor- 
phisms of the unit circle modulo the rotations |KMj . etc. 

In contrast with all the examples we gave above, the energy of the mod- 
ified Euler-Poisson equation (mEP) is not a quadratic form, therefore it can 
not be interpreted as a Riemannian metric. However, there still is a vari- 
ational problem on the cotangent space of the configuration space of this 
equation. Here we derive the equation (mEP) from this variational prob- 
lem. Note that all the computations that follow are formal. 

Let g be a Lie algebra with the bracket operation [•, •] and g* be its dual 
given by the pairing 

(•,•) :0* xg^M. 



Then g* with the Lie-Poisson bracket defined by 



{F,G}{m) = {m, 



6F_ 5G_ 
5m ' 5m 



) 



for G : g* — > M is a Poisson manifold. 

Hamilton's equations for H € C°°{P) on a Poisson manifold (P, {•,•}) 
are given by 

F = {F, H] for ah F G C7°°. (26) 
Given a Lie-Poisson manifold, Hamilton's equations can be written as 



F{m) = {F,H}{m) 

5H 6F 



{^,dtm) = -{m 
dm 



5m 5m 



Therefore 



dtm = -adljj/sm^ (27) 

is an equivalent formulation of the Hamilton's equation on a Lie-Poisson 
manifold. 
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Here we exploit the tools and techniques used to study the Hamilto- 
nian formulation of the Euler equations for a compressible fluid [MRW ^ to 
show that the modified Euler- Poisson equation (mEP) can be derived as a 
Hamiltonian equation. 

On the Cartesian product space Diff(T™) x C°°(T'") of the group of 
diffeomorphisms of T™ and the vector space C°°(T'") of all smooth functions 
on T™', the operation 

((/>, a) o (^/;, 6) = (0 o ^/>, a o + b), 

called the semidirect product induces a Lie group structure. We denote this 
group by 

G = Diff(T™) X C°°(T") 

following the conventional notation for semidirect product spaces. The cor- 
responding Lie algebra is the space q = Vect(T™) x C°°(T™') with the bracket 

[{v,a),{w,b)] = {[v,w],Cyb- C^a) 

where v,w (z Vect(T'") and a,b ^ C°°(T'"). Here [v,w] is the usual com- 
mutator of vector fields on T™ and CwCi is the Lie derivative of a in the 
direction of w and is given by >C^a = ^j^^Q (ao^*) where is any curve on 
Diff(T'^) such that C*|s=o = and ^|g_Q C'^ = w. Note that in this setting 
the composition of the diffeomorphisms is the group operation on Diff (T™"), 
the composition of a smooth function with a diffeomorphism a o is the 
natural action of the diffeomorphism 7 on the function a. In general, the 
semidirect product structure on the Cartesian product of a Lie group and a 
vector space on which the group acts is defined using the group operation 
and the action of the group on the vector space (see |AKj . |MRW _ ) . 

In this context, a Hamiltonian formulation of the modified Euler- Poisson 
equation (mEP) can be stated as follows: 

Theorem 4. The modified Euler-Poisson equation (mEP) is a Hamiltonian 
equation on q* with respect to the linear Lie-Poisson structure and the energy 
function 

H{M,n) = - ( — {M,M) +^{n)] dx (28) 
Jt^ \2n J 

where M = nv e Vect(T™-) and ^'{n) = A-'^{n). 

Proof. We want to derive the equations for v and n from 

dtm = -adlH/s^m (29) 
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where m = (M,n) = (nv,n) and {v,n) € Q* . 

The variational derivative is given by ( |^ , j with 



6Hi 

Jm 



-M/n = -V, 



(30) 
(31) 



^ = ^(M,M) - cl>'(n) = ^{v,v) - <I>'(n). 
Evaluating equation ()29() on an arbitrary pair (if, 6) € 0, we obtain 

{dtm, {w, b)) = -{adljj^i^^m, {w, h)). 

Then by the definition of the coadjoint operator ad* and the bracket on g, 
we have 



{dtm,{w,b)) = -{iM,n) 
= ((M,n) 



w 



6Hi 6Hi 
5Hi 



{w,b) 



) 



r r h\\ 

Lw— LsHiO ]). 

On TW J 



In what follows we identify the dual space q* with the algebra g using the 
pairing (•, •) on g* x g given by 



(32) 



{{v,a), {w,b)) = / {v,w) dx + ab dx. 

Jjm Jjm 

Then substituting nv for M and using pOj) and ()31() . we obtain 

{dtm,{w,b)) = / {[w, —v],nv)dx+ 

+ n{c^{]^{v,v) -^'{n))- C^b^dx. 
By the definition of the bracket [•, •] on Vect(T™), we have 
{dtm,{w,b)) = / {{v.V)'w,nv) — {{w.V)v,nv)dx 

+ n(^Cyj{^{v,v) -<P'{n))-£-yb^dx. 
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Furthermore, we can compute the Lie derivatives on the right hand side and 
write the above equahty as fohows 

(dtm, {w,b)) = / {{v.V)'w,nv) — {{w.V)v,nv)dx+ 

Jfm 

+ I {nw,grad{-{v,v) — <I>'(n))) + {nv,giad{b))dx. 

Jfm 2 

Using the identities 

{v, grad{v, w)) = {w, {v.'V)v) + {v, {v.'V)w) 

and ^ 

{w, -grad(z;, v)) = {v, {w.V)v) 

we obtain 

{dtm,{'w,b)) = / {giad{v,w),nv) — {{v.V)v,nw) — 
JJm 

— (nw, grad($'(n))) + (nw, grad(6))dx. 

Integrating by parts the first and the last summands on the right hand side, 
we obtain 



{dtm, {w,b)) = I —diY{nv){v,w) — {{v.V)v,nw) — {nw,grad{^'{n)))—bdiY{nv)dx. 
Note that 



{dtm, {w, b)) = {{{dtn)v + n{dtv),dtn), {w, b)). 
Then by ((SH), we have 

dtn = — div(nt;), 
{dtn)v + n{dtv) = —dvv{nv)v — n{v.V)v — ngrad($'(n)) 



which is equivalent to 



dtn = — div(nz;), 

dtv = —{v.\/)v — grad(<I>'(n)). 



□ 



Note that for one space dimension (m = 1) the hamiltonian Hi in ()28|) 
is given by 

Hi = j]^{v^n + {K-^d^nf + {K-^nf)dx (33) 
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in terms of v and n. Then using the differential operator Vi that is defined 
as 

one can rewrite the equation (mEP) in the hamiltonian form 

The Poisson bracket induced by the matrix differential operator Vi is given 
by 



{F,H}{v,n) = {{v,n) 



' (5F 


5F\ 


fSH 


6H\' 




Sn J 




6n J 



) 



SF^ 6H 5H ^ 6F\ , 

-^Ox-r ^(^x^ I ax. 

on dv on ov 



Another conserved quantity for (mEP) is 

H2 = nv dx. (35) 



For one space dimension (m = 1) we can use H2 to write (mEP) in yet 
another form as 

where V2 is defined as 



2 



dxV —(nd + dn) 



We prove theorem |21 by showing that H34|) and (jSHJ are Hamiltonian forms 
of the modified Euler-Poisson equation (mEP) and that the induced Poisson 
structures are compatible. 

proof of theorem |3J The matrix differential operator Vi is skew- 
adjoint and does not depend on v nor n nor any of their derivatives, therefore 
the bracket given by Pi satisfies the Jacobi identity hence is indeed a Poisson 
bracket. 

We can easily check that D2 is skew-adjoint as well: 

h,(l)2)V2{9i,e2)dx 
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01 (A ^9^6*1 - 92dj:v) + 4>2{0idxV - ndx02 - dx{ne2))dx 

To verify the Jacobi identity for the bracket induced by P2 we adapt the 
notation of prolongations (see for details). Let 62 be the functional 
bivector associated to V2: 

Q2 = l Jiei,92)V2i9i,e2)dx 

= ^ y {^1 A A-^dx6i -01 A e2dxv 
+02 A Oidr^v - 6*2 A ndx92 - O2 A dx{92n)}dx 
= \ /{^i ^ ^'"^dJi - 2{dxv)9i A 92 - 202 A ndj2}dx. 
Then P2 is Hamiltonian since 

pr -^vA^2) = -\j {2{9Av - ndj2 - 9.(^2^)) A ^2 A 8^2 

+da^{A-^dx9i - 92dxv) A A 92]dx 
= - j {Pr^dl9i A 01 A 92]dx 

= - j {A-26)i A 01 A 92]dx 

= -\j {dli^^'^Oi A A-'^9i) A 92}dx 
= 0. 

These two Hamiltonian structures, (jHlJ and (jHSJ, are compatible, i.e. 
the equation (mEP) is bihamiltonian in one space dimension. To prove the 
compatibility it is enough to check that 

pr vxi^eiBvi) + pr ^VidCdvi) = (37) 
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holds where Q-jji denotes the corresponding bivector for Pj . Both summands 
in H37() vanish: 

pr ^rj,,e{&v,) = \ j {-2dJi A 02 A 8,02 - 28% A0iA 02}dx 

= - j d,{0,A02A8j2)dx 
= 0, 

and similarly we have 

pr v-p^eiOvi) = 

since 

®v^=^J {-Oi A 8,02 -O2A d,0i}dx 

= -J{0^A 8j2}dx. 

Therefore the modified Euler-Poisson equation (mEP) is bihamiltonian for 
m = 1 and this completes the proof of theorem |51 

□ 
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